Calculation of biological effectiveness of SOBP proton beams: a TOPAS Monte Carlo study.

Objective.This study aims to investigate the biological effectiveness of Spread-Out Bragg-Peak (SOBP) proton beams with initial kinetic energies 50-250 MeV at different depths in water using TOPAS Monte Carlo code.Approach.The study modelled SOBP proton beams using TOPAS time feature. Various LET-based models and Repair-Misrepair-Fixation model were employed to calculate Relative Biological Effectiveness (RBE) for V79 cell lines at different on-axis depths based on TOPAS. Microdosimetric Kinetic Model and biological weighting function-based models, which utilize microdosimetric distributions, were also used to estimate the RBE. A phase-space-based method was adopted for calculating microdosimetric distributions.Main results.The trend of variation of RBE with depth is similar in all the RBE models, but the absolute RBE values vary based on the calculation models. RBE sharply increases at the distal edge of SOBP proton beams. In the entrance region of all the proton beams, RBE values at 4 Gy i.e. RBE(4 Gy) resulting from different models are in the range of 1.04-1.07, comparable to clinically used generic RBE of 1.1. Moving from the proximal to distal end of the SOBP, RBE(4 Gy) is in the range of 1.15-1.33, 1.13-1.21, 1.11-1.17, 1.13-1.18 and 1.17-1.21, respectively for 50, 100, 150, 200 and 250 MeV SOBP beams, whereas at the distal dose fall-off region, these values are 1.68, 1.53, 1.44, 1.42 and 1.40, respectively.Significance.The study emphasises application of depth-, dose- and energy- dependent RBE values in clinical application of proton beams.

Microdosimetry-based investigation of biological effectiveness of252Cf brachytherapy source: TOPAS Monte Carlo study.

Objective.To investigate biological effectiveness of252Cf brachytherapy source using Monte Carlo-calculated microdosimetric distributions.Approach.252Cf source capsule was placed at the center of the spherical water phantom and phase-space data were scored as a function of radial distance in water (R= 1-5 cm) using TOPAS Monte Carlo code. The phase-space data were used to calculate microdosimetric distributions at 1µm site size. Using these distributions, Relative Biological Effectiveness (RBE), mean quality factor (Q̅) and Oxygen Enhancement Ratio (OER) were calculated as a function ofR.Main results.The overall shapes of the microdosimetric distributions are comparable at all the radial distances in water. However, slight variation in the bin-wise yield is observed withR. RBE,Q̅and OER are insensitive to R over the range 1-5 cm. Microdosimetric kinetic model based RBE values are about 2.3 and 2.8 for HSG tumour cells and V79 cells, respectively, whereas biological weighting function-based RBE is about 2.8. ICRP 60 and ICRU 40 recommendation-basedQ̅values are about 14.5 and 16, respectively. Theory of dual radiation action based RBE is 11.4. The calculated value of OER is 1.6.Significance.This study demonstrates the relative insensitivity of RBE,Q̅and OER radially away from the252Cf source along the distances of 1-5 cm in water.

Comparison of 126 MeV antiproton and proton-a FLUKA-based microdosimetric approach.

Objective. This study aims at comparing dosimetric parameters of 126 MeV antiprotons and protons using microdosimetric approach.Approach. Microdosimetric distributions of 126 MeV proton and antiproton beams at 1µm site size are calculated using the Monte Carlo-based FLUKA code. The distributions are calculated at various depths along the central axis in water phantom as well as at different off-axis locations. The study also includes calculations of secondary radiations produced by antiprotons and protons. Mean quality factor,Q¯is calculated using the ICRP 60 and ICRU 40 recommendations. The Relative Biological Effectiveness (RBE) of HSG tumour cell at 10% survival level is calculated based on Microdosimetric Kinetic Model.Main results.Q¯ICRP,Q¯ICRUand RBE for antiprotons are higher by a factor of about 3.60, 3.41 and 1.24, respectively, at Bragg-peak and higher by a factor of about 1.41, 1.76 and 1.05, respectively, at plateau region of depth-dose profile when compared to protons. At 15 cm depth along the central axis,Q¯ICRP,Q¯ICRUand RBE for protons are higher by a factor of about 1.42, 1.66 and 1.26, respectively, when compared to antiprotons. At the off-axis distance (Ld) of 6 cm (at 11.5 cm depth in water),Q¯ICRPandQ¯ICRUof protons are higher than that of antiproton whereas the trend is opposite at off-axis distance of 4 cm. AtLd= 4 cm (at 11.5 cm depth in water), RBE of antiprotons is higher by about 4% than protons whereas atLd= 6 cm, RBE of protons is higher by about 13% than antiprotons.Significance. The study shows that antiproton radiotherapy is advantageous as compared to protons considering enhancements in the absorbed dose and RBE-weighed dose values at the Bragg-peak.

Surface dose rate variations in planar and curved geometries of 106Ru/106Rh plaque sources for ocular tumors.

PURPOSE: Investigation of surface dose rate variation with respect to the source configuration of 106Ru/106Rh eye plaque. To explore an alternate way to determine activity of brachytherapy plaques. METHODS: The surface dose rates of 106Ru/106Rh plaque developed indigenously were measured by extrapolation chamber. To rule out possibility of any error in the activity distribution and quantity, same source was used in two different configurations namely planar and curved. EBT3 Gafchromic film was used for determination of uniformity in activity. Monte Carlo-based Codes EGSnrc and FLUKA were used to calculate dose rate in tissue, percentage depth dose and for determination of activity. Parameters and correction factors were estimated using simulations. RESULTS: The measured reference absorbed dose rates for planar and curved 106Ru/106Rh eye plaques are found to be 589 ± 29 mGy/h and 560 ± 28 mGy/h, respectively. The difference in the reference absorbed dose rate of curved eye plaque is about ~5% as compared to planar configuration. The FLUKA-calculated dose values are almost independent of cavity length of the extrapolation chamber for both eye plaques. The FLUKA-based dose rates per µCi 106Ru/106Rh are about 17.28 ± 0.08 mGy/h and 16.48 ± 0.06 mGy/h, respectively for planar and curved eye plaques which match well with the measurements. The calculated activities for planar and curved eye plaques are 34.08 µCi and 33.98 µCi, respectively. CONCLUSIONS: Surface dose rates for a prototype 106Ru/106Rh eye plaque with different configurations were estimated using simulations and measured experimentally. An alternate way to determine activity of beta-gamma brachytherapy plaque has been proposed.

Microdosimetry-based relative biological effectiveness calculations for radiotherapeutic electron beams: a FLUKA-based study.

Based on FLUKA, the present study is aimed at calculating the microdosimetric distributions of electron beams (6, 12 and 18 MeV) for radiotherapy as a function of depth in water at a site-size of 1 µm using a tissue-equivalent proportional counter (TEPC). Using the calculated microdosimetric distributions, the depth-specific relative biological effectiveness (RBE) of electron beams used in radiotherapy is calculated based on the theory of dual radiation action (TDRA) and the microdosimetric kinetic model (MKM). The TDRA-based calculation shows the variation of RBE of an electron beam with the absorbed dose and depth in water. In this study, we compared the RBE values calculated based on the TDRA and MKM. The FLUKA-based microdosimetric distributions in water obtained using the pre-calculated electron fluence spectra resulted in an improvement in the computational efficiency by a factor of 110 when compared with a full simulation. Depending on the beam energy and depth of water, RBETDRA was in the range 0.67-0.78. RBEMKM at 10% survival of HSG tumor cells was 0.84, which was nearly independent of the depth and beam energy.

TISSUE-EQUIVALENCE OF H2 GAS FOR MICRODOSIMETRY IN NEUTRON FIELDS: A GEANT4 MONTE CARLO STUDY.

Hydrogen (H2) as filling gas in Tissue-equivalent proportional counter (TEPC) for measurements of microdosimetric distributions in neutron fields is investigated using the Monte Carlo-based Geant4 toolkit. The neutron fields considered are monoenergetic neutrons (1 keV-14 MeV) and ISO reference neutron sources 241Am-Be, 241Am-B, 252Cf and 252Cf + D2O. Based on the energy deposited in the gas cavities (tissue-equivalent (TE) propane and H2) of the TEPC, density scaling correction, ${f}_{\rho }$, to be applied on H2 gas to achieve tissue-equivalence is calculated for the above sources at 1µm site size. The calculated value of ${f}_{\rho }$ for the ISO-neutron sources is 0.40, which suggests that energy deposition in the gas cavities is predominantly due to crossers. In the case of monoenergetic neutrons, depending upon the energy, ${f}_{\rho }$ is in the range of 0.11-0.45. The study shows that the TE propane and H2-based microdosimetric distributions are comparable when the density of H2 is scaled appropriately. Mean quality factor $\overline{Q}$ calculated based on microdosimetric distribution increases initially with neutron energy up to 100 keV and thereafter decreases with energy. Depending upon the neutron energy, TE propane and H2-based $\overline{Q}$ values are comparable to within 3-15%, and for the ISO sources, the comparison is within 5-8%. For the ISO neutron sources, $\overline{Q}$ values are in the range of ~10-15.

APPLICABILITY OF PURE PROPANE GAS FOR MICRODOSIMETRY AT BRACHYTHERAPY ENERGIES: A FLUKA STUDY.

Applicability of pure propane gas for microdosimetric measurements at photon energies relevant in brachytherapy is studied using the Monte Carlo-based FLUKA code. Monoenergetic photons in the energy range of 20-1250 keV and brachytherapy sources such as 103Pd, 125I, 169Yb, 192Ir, 137Cs and 60Co are considered in the study. Using the calculated values of energy deposited in the sensitive region of LET-1/2 tissue-equivalent proportional counter filled with pure propane gas and tissue-equivalent propane gas, values of density scaling factor for the site sizes of 1 and 8 µm are obtained. The study shows that density of propane gas should be lowered by a factor of about 0.93 for 169Yb, 192Ir, 137Cs and 60Co sources for the site sizes of 1-8 µm. For 125I source, the density of propane gas requires a scaling of 0.93 for 1 µm site size, whereas for site sizes 2-8 µm, density need not be altered. 103Pd source does not require density scaling for site sizes 1-8 µm.

MONTE CARLO-BASED INVESTIGATION OF MICRODOSIMETRIC DISTRIBUTION OF HIGH ENERGY BRACHYTHERAPY SOURCES.

FLUKA-based Monte Carlo calculations were carried out to study microdosimetric distributions in air and in water for encapsulated high energy brachytherapy sources (60Co, 137Cs, 192Ir and 169Yb) by simulating a Tissue Equivalent Proportional Counter (Model LET1/2) having sensitive diameter of 1. 27 cm for a site size of 1 µm. The study also included microdosimetric distributions of bare sources. When the sources are in air, for a given source, the source geometry does not affect the y¯F and y¯D values significantly. When the encapsulated 192Ir, 137Cs and 60Co sources are in water, y¯F and y¯D values increase with distance in water which is due to degradation in the energy of photons. Using the calculated values of y¯D, relative biological effectiveness (RBE) was obtained for the investigated sources. When 60Co, 137Cs and 192Ir sources are in water, RBE increases from 1.03 ± 0.01 to 1.17 ± 0.01, 1.24 ± 0.01 to 1.46 ± 0.02 and 1.50 ± 0.01 to 1.75 ± 0.03, respectively, when the distance was increased from 3-15 cm, whereas for 169Yb, RBE is about 2, independent of distance in water.

INVESTIGATION OF APPLICABILITY OF PURE PROPANE GAS FOR MICRODOSIMETRY AT NEUTRON FIELDS: A MONTE CARLO STUDY.

Applicability of pure propane gas for microdosimetric measurements in neutron fields was investigated using the FLUKA Monte Carlo code. Monoenergetic neutrons in the energy range 1 keV-20 MeV and the ISO-neutron sources such as 241Am-Be, 241Am-B, 252Cf and 252Cf + D2O were considered in the present study. The tissue-equivalent proportional counter (TEPC) simulated in the study was LET-1/2 (by Far West Technology) with site sizes 1, 2 and 8 µm. The study demonstrates that for a given site size, the TEPC filled with tissue-equivalent propane and pure propane gases produce similar microdosimetric distributions when the density of pure propane gas is lowered appropriately. For the ISO-neutron sources, the density of propane gas requires scaling by a factor 0.85. For the monoenergetic neutrons, depending upon the neutron energy, the values of scaling factors are in the range of 0.58-0.93.

FLUKA-BASED MONTE CARLO INVESTIGATION OF MICRODOSIMETRIC DISTRIBUTIONS OF TELECOBALT BEAM.

FLUKA-based Monte Carlo calculations of microdosimetric distributions in water phantom involving a walled spherical Tissue-Equivalent Proportional Counter filled with tissue-equivalent propane gas have been studied for an indigenously developed telecobalt machine. The simulated site size considered in the study was 2 µm. In the Monte Carlo calculations, field size was varied from 10 cm × 10 cm to 35 cm × 35 cm and the depth was varied as 5-20 cm. The study also includes calculation of microdosimetric distributions with a 30° wedge filter. The efficiency of the calculations was improved up to a factor of 26 by choosing appropriate cut off values for production and transport of electrons. The calculated microdosimetric distributions of telecobalt machine is distinctly different from that of a bare 60Co source which is attributed to the influence of scattered photons from the machine head and the water phantom.